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ABSTRACT 

Swift J1357.2-0933 is one of the shortest orbital period black hole X-ray transients (BHTs). 
It exhibited deep optical dips together with an extremely broad Ha line during outburst. We 
present 10.4-m GTC time-resolved spectroscopy during quiescence searching for donor star 
absorption features. The large contribution of the accretion flow to the total luminosity pre¬ 
vents the direct detection of the companion. Nevertheless, we constrain the non-stellar con¬ 
tribution to be larger than ~ 80% of the total optical light, which sets new lower limits to the 
distance (d > 2.29kpc) and the height over the Galactic plane (z > 1.75kpc). This places the 
system in the galactic thick disc. We measure a modulation in the centroid of the Ha line 
with a period of P = 0.11 ± 0.04d which, combined with the recently presented FWHM-X 2 
correlation, results in a massive black hole (Mi > 9.3 Mq) and a ~ M2V companion star 
(M 2 0.4M©). We also present further evidence supporting a very high orbital inclination 

(1^80°). 
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1 INTRODUCTION 

Galactic X-ray binaries offer the best opportunity to detect stel¬ 
lar mass black holes (BHs) through the accretion luminosity pro¬ 
duced by material from the companion star falling into its grav¬ 
itational well. The vast majority of BHs are detected in transient 
low-mass X-ray binaries (LMXBs), which spend most part of their 
lives in a faint, quiescent state (about lO^^^^^ergs^^; see e.g. 
Armas Padilla et al. 2014b). They are discovered during occasional 
outburst, where their X-ray luminosity increases above ~ 10 per 
cent of the Eddington limit. There are only 17 dynamically con¬ 
firmed BHTs, as well as ~ 33 candidates (see Casares & Jonker 
2014; Corral-Santana et al. 2015). 

Swift J1357.2-0933 (hereafter J1357) is a LMXB X-ray tran¬ 
sient with galactic coordinates I = 328.702° and b = 4-50.004° 
(Krimm et al. 2011), whose peak X-ray luminosity place it in the 
very faint regime (see Armas Padilla et al. 2013 for an X-ray anal¬ 
ysis during outburst). The detection of optical dips both in out¬ 
burst (Corral-Santana et al. 2013, hereafter CS13) and quiescence 
(Shahbaz et al. 2013) suggests J1357 is a very high inclination sys¬ 
tem. The orbital period is among the shortest of its class (P = 2.8h, 
CS13) and the mass function is constrained to be /(Mi) > 3 .OM 0 , 
where the radial velocity of the donor (fe) was estimated from 
the double-peak separation of the Ha profile (Orosz et al. 1994; 
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Orosz & Bailyn 1995). This strongly advocates for the presence 
of a BH. Rau, Greiner & Filgas (2011) proposed a tentative dis¬ 
tance to the system of 4 ~ 1.5kpc considering the donor star to 
be the dominant source of quiescent light in the optical and near- 
infrared regime. However, the quiescent spectral energy distribu¬ 
tion (SED) is best described by a single power-law model, which 
suggests little, if any, thermal contribution from the secondary star 
(Shahbaz et al. 2013). This implies that only lower limits to the dis¬ 
tance can be obtained unless the donor star contribution is properly 
characterized. 


2 OBSERVATIONS 

J1357 was observed with the Optical System for Imaging and 
low-Intermediate-Resolution Integrated Spectroscopy (OSIRIS) lo¬ 
cated in the Nasmyth-B focus of the 10.4-m GTC, in La Palma 
(Spain). We used the R500R optical grism combined with a 1.0” slit 
(R = 352, 4.88 A/pix), covering the spectral range 4800— 10000 A. 
Observations were obtained in four different nights within April to 
June 2014. 

Eleven consecutive spectra per night were acquired at air- 
masses in the range ~ 1.3 — 1.7, with individual exposure times of 
875s (~ 3h per night) sampling the proposed 2.8h orbital period. 
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Figure 1. Averaged spectnjm of J1357. Accretion disc emission features (Ha and He I) are marked with blue, dotted vertical lines. Regions contaminated by 
sky subtraction residuals are plotted as red, doted vertical lines. The expected M-type donor star absorption features are shown as green bands (TiO molecular 
bands) and green, dotted vertical lines (Nal and Call triplet). 


We used IRAF' standard routines for bias and flat-field corrections. 
The slit was rotated to PA = —53.16° in order to allow for simulta¬ 
neous observations of the object and a nearby early-type dwarf star 
located at 2.06arcmin NE of the target. The wavelength calibration 
was obtained from observations of calibration arcs and instrumen¬ 
tal flexure was calculated from the drift of the [O I] 5577 A sky line 
and subsequently corrected in every individual spectrum. 

The spectrum of the calibration star is classified as a F8V 
after comparison with templates from the MILES spectral library 
(Sanchez-Blazquez et al. 2006, Falcon-Barroso et al. 2011) using 
minimization routines within MOLLY (e.g. Casares et al. 1996). 
In order to maximize the chances of detecting donor star absorption 
features in J1357, we produced a set of telluric spectra by remov¬ 
ing the F8V spectral lines, which were subsequently substracted 
from the target to produce spectra free of telluric features (e.g. 
Beekman et al. 2000). 

The above strategy yielded four sets of eleven spectra. Points 
deviating more than 3(7 above or below the continuum, probably 
caused by bad pixels or residuals from the sky subtraction, were 
interpolated. One spectrum was discarded on each of the two first 
nights due to very poor signal-to-noise ratio (S/N) caused by the 
presence of clouds. Absorption features from the donor star are not 
evident in any of the spectra, nor in the combined spectrum (see 
Fig. 1). 


3 RESULTS 

Using the 42 GTC spectra we have determined the veiling factor, 
the orbital period and the main properties of the Ha line in quies¬ 
cence. 

3.1 Search for companion star features: Skew-mapping 

The skew-mapping technique has been applied in cata¬ 
clysmic variables to detect weak companion star features (e.g. 
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Vande Putte et al. 2003, Smith et al. 2005). Our analysis compares 
a template M3V star spectrum with different Doppler-corrected 
averaged spectra created from our database. Each averaged 
spectrum is produced by considering a specific zero phase in the 
range 0-1 (we inspect all values with a 0.05 step) and velocity of 
the companion star between K 2 = 0 —2000kms^*, in 100kms^* 
steps. Individual spectra are shifted and coadded assuming the 
orbital period of CS13. The cross-correlation with the M3 template 
star is expected to reveal a significant peak above noise at the 
correct parameters for the companion star motion. This technique 
did not reveal any preferred value in the A '2 - zero phase parameter 
space, not even after considering several different periods around 
the one proposed in CS13 (0.09 — 0.15d^^, steps of O.OOld^*). 
Therefore, we conclude that no companion star features are present 
in our spectra. 


3.2 Veiling factor 

Magnitudes in the OSIRIS / filter were obtained from the acqui¬ 
sition images of each night: /[ = 21.71 ±0.08; (2 = 21.45 ± 

0.06; r^ = 20.84±0.08; = 21.00±0.11. The first two values 

(April 29th and 30th) are similar to those reported in Shahbaz et al. 
(2013, r' = 21.54 ±0.35). On the other hand, the latest two nights 
(June 2nd and 28th) reveal a somehow brighter system (but still 
in quiescence). We note that although this might be caused by 
the known strong short-term variability displayed by the system 
(Shahbaz et al. 2013), the spectra of the first two nights are clearly 
noisier under similar sky conditions, reflecting a true drop in bright¬ 
ness. 

The absence of companion star features imposes a minimum 
constraint to the veiling factor (Y), defined as the fractional con¬ 
tribution of the accretion related luminosity (La^c) 1° l^e total flux 
(Laac + ^2 ’ where ±2 is the donor star luminosity) in the OSIRIS 
r' filter wavelength range (5500 — 7400 A). On the other hand, by 
considering the relation between the orbital period and the mean 
density presented in Faulkner, Flannery & Warner (1972), and tab¬ 
ulated values for main-sequence stars (Cox 2000), the spectral type 
of the companion star is constrained to be M2V or later. We com¬ 
pared our spectrum with templates of main-sequence stars from the 
MILES Spectral Library. For each spectral type, we measured the 
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normalized flux of the deepest photospheric absorption line present 
in the template spectrum. The veiling necessary to make these fea¬ 
tures shallower than the noise level (3(7) within the corresponding 
spectral region of the J1357 spectrum is a lower limit to the veiling 
factor (see Mata Sanchez et al. 2015 for a similar analysis of NIR 
spectra). 

We find that the required veiling factor for a M2V star to be 
swamped by the accretion flux isX >0.81. This limit is also valid if 
the donor is evolved and its spectral type is later than M2V (we have 
explored spectral types as late as M6V ). This result is obtained 
using the averaged quiescent spectrum of the first two nights where 
the system is fainter, and therefore, a higher contribution of the 
companion star is expected. 

We have explored spectral types as late as M6V, obtaining 
similar results than for M2V. We note that using slightly different 
tabulated values for dwarf star mean densities (e.g. Carroll & Ostlie 
2006), does not affect our results. 

3.3 Orbital period from a two-Gaussian model fit 

We fit the Hot double-peaked emission profile with a model consist¬ 
ing of two Gaussians with equal FWHM. The height of the Gaus- 
sians and offset with respect to the Hot rest wavelength (6562.78 
A) are left as free parameters. 

Ha centroid velocity 

The mean of the Gaussians offsets is expected to trace the motion of 
the accretion disc around the centre of mass of the system. We note 
that the mean offset value for each night was significantly different 
(see Fig. 2b). This behaviour has been reported in other BHTs, for 
example in XTE Ill 18-1-480 it is explained as a result of disc pre¬ 
cession (Torres et al. 2002). Recently, Torres et al. (2015) detected 
variations in the systemic velocity of J1357. They did not consider 
this result to support the precessing disc scenario, but interpreted 
as systematic effects. However, our longer database suggests that 
instead this is a consequence of a physical process taking place in 
the system. In order to fit the whole sample, we then decided to 
de-trend the mean nightly offsets. 

We performed a Lomb-Scargle normalised periodogram of 
the (de-trended) mean offset variations and the result is pre¬ 
sented in Fig. 3. We used the period time-series analysis package 
(Dhillon, Privett & Duffey 2001). A Gaussian fit to the broad peak 
with the highest power yields a period of P = 0.11 ± 0.04 d (1 a). 
This is consistent with the one measured by CS13 using a longer 
database with higher S/N during outburst (P = 0.117 ± 0.013d). 
We find two families of peaks consistent with the outburst period. 
Gaussian fits result in P = 0.109 ± 0.005 d and P = 0.121 ± 0.006d, 
respectively. The former is favoured by its highest amplitude (de¬ 
picted as a red, stripped band in Fig. 3). Nevertheless, further obser¬ 
vations are necessary to confirm this result. We note that compara¬ 
ble results are obtained by performing sinusoidal fits and minimise 
chi-square. 

The data were folded on the previously obtained period and 
are shown in Fig. 2b, where different colours and symbols are used 
to mark different systemic velocities observed on each night. We 
subsequently performed a non-linear least squares fit to the data 
using the following sinusoidal function: 

V = y+Ki sm(2jt{(j> - (j> o)) 

We obtained Ki = 44 ± 3kms^*, y = 75 ± 2kms^* and (/> o = 



Phase 


Figure 2. Orbital evolution of several Hc« parameters obtained from our 
two-Gaussian model fit (three top panels) and from a diagnostic diagram 
(bottom panel) folded with the best fit period obtained from periodogram 
analysis. Panel a: intensity ratio between red and blue peak. Best fit is rep¬ 
resented as a black, solid line. Panel b: radial velocities of the non-detrended 
centroid as measured by the two-Gaussian model. In both panels, each night 
is represented by different symbols and colours as follows: red circles (April 
29th), blue triangles (April 30th), green upside-down triangles (June 2nd) 
and orange squares (June 28th). The best sinusoidal fit to the non-detrended 
data is also shown. Panel c: radial velocities of the line centroid obtained 
from 10 phase bin averaged spectra after de-trending. The best sinusoidal 
fit is shown as a blue, solid line. Panel d: same as previous panel but using 
radial velocity curves from the diagnostic diagram (Gaussian separation of 
a = 4500kms^'). Data are plotted over two phase cycles for clarity. 

0.71 ±0.01. Even if it might be qualitatively acceptable, the high 
reduced chi-square value (;); ^=17.26) shows that the non de-trended 
data does not properly fit to a sinusoid, and de-trending is nec¬ 
essary for an optimal analysis. Subsequently, the original spectra 
were combined into 10 phase bins after removing each night’s mean 
value. The fit results in o = 0.720±0.009 and K\ = 61 ±4kms^* 
for (see Eig. 2c). 

Intensity ratio and double peak separation 

The intensity ratio of the peaks of the Htz profile (defined as 
/jed/7blue) exhibits a periodicity similar to that observed in the cen¬ 
troid offsets. Folding these data onto the previously obtained period 
and fitting a sinusoid results in ^=5.12 (see Fig. 2a). Torres et al. 
(2015) reported asymmetric Ha profiles, but no recurrent vari- 
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to standard deviations). This result is consistent with previous de¬ 
terminations of the systemic velocity in quiescence (y= —130± 
50kms^*; Torres etal. 2015) and the compact object’s velocity 
reported in outburst {K\ = 43±2kms^*; CS13). However, the 
systemic velocity value obtained from the outburst periodogram 
( 7 ~ —150kms^*, CS13) does not seem to be in agreement with 
our results. 

We note that, in contrast to the two-Gaussian profile modelling 
(see Section 3.3), nigthly de-trending is not necessary for the radial 
velocity curves obtained with the double Gaussian technique. This 
can be explained by considering that the two-Gaussian modelling 
traces the core of the Ha double peak profile, and therefore, outer 
parts of the disc where disc preccesion effects might be important. 
However, the diagnostic diagram is sensitive to the emission line 
wings, which traces inner parts of the disc less affected by com¬ 
panion star tidal forces. 


Figure 3. Periodogram obtained with the Lomb-Scargle method. Red, solid 
line shows a Gaussian fit to the highest group of peaks. The red, dot¬ 
ted horizontal line represents the Iff uncertainty (P = 0.11 ±0.04d). We 
inspected frequencies up to 60 d“', but the power spectrum for frequen¬ 
cies over 30 d“* is nearly fiat. The red, dashed vertical line refers to the 
strongest peak consistent with the outburst period proposed by CS13. Grey 
band depicts the constraint to the period (P = 0.117 ±0.013d) from out¬ 
burst data (CS13). The red, stripped band refers to the standard deviation 
obtained from a Gaussian fit to the family of peaks inside the grey band 
(P = 0.109 ± 0.005 d). 

ability was observed. The detection of this periodic evolution of 
the profile along four different nights suggest an scenario with an 
asymmetric outer disc structure, perhaps a hot-spot on a tidal arm. 
This would also explain the deviation observed from a perfect fit, 
since other processes affecting the intensity ratio might be at play. 

The double-peak separation does not show significant variabil¬ 
ity, and it is consistent with a constant value. The weighted average 
of the double-peak separation is Dp = 2430 ± 50kms^*, signifi¬ 
cantly larger than the one measured in outburst by CS13 {Dp = 
1790±67kms^*). This implies that the outer disc velocity (v^ = 
1215±25kms^*) is higher in quiescence, which is an expected be¬ 
haviour since the disc expands during outburst, reaching areas with 
lower velocities. The reported value is consistent within 2(7 with 
that reported by Torres et al. (2015, Dp = 2340±20kms^*). 

3.4 The diagnostic diagram 

The broad wings of the Ha profile are supposed to trace the motion 
of the compact object. We apply the diagnostic diagram (aka double 
Gaussian technique Shafter, Szkody & Thorstensen 1986), aiming 
at constraining radial velocity variations of the Ha wings. We in¬ 
spected several Gaussian separations {a = 2500-- 5000kms^*) to 
perform a diagnostic diagram analogous to that presented in CS13. 
Unfortunately, our limited database and the lower S/N ratio pre¬ 
vents us from obtaining conclusive determinations for the A"! and 
Y values. A wide plateau on the evolution of the fitted parameters 
from a = 4000kms^*, combined with the inspection of the indi¬ 
vidual radial velocity curves, suggests a = 4500kms^^ as the pre¬ 
ferred value (see the folded, radial velocity curve for this separation 
in Fig. 2d). 

Mean values of the compact object orbital velocity {Ki = 
40 ± 12kms^*) and systematic velocity of the system (y = —79 ± 
13kms^*) are obtained from sinusoidal fits to radial velocity 
curves in the range a = 4000 — 5000kms^* (uncertainties refer 


3.5 H-alpha full-width-at-half-maximum 

We have measured the full-width-at-half-maximum (FWHM) of 
the Ha emission profile in each individual spectrum. After sub- 
stracting quadratically the instrumental resolution, the average 
value is FWHM = 4152 ± 209 km s^*, where the uncertainty re¬ 
flects the standard deviation of the 42 measures. This value is larger 
than the one measured in outburst (FWHM ~ 3300kms^*, CS13) 
and consistent with the quiescence value reported in Torres et al. 
(2015) and Casares (2015), i.e. FWHM = 4025 ± 1 lOkms^' and 
FWHM = 4085 ± 328 km s^ * respectively. The previously reported 
broadest FWHM corresponds to XTE HI 18-1-480 (FWHM ~ 
2500km s^ *, Torres et al. 2004). 


4 DISCUSSION 

Observations taken during the outburst decay have sometimes re¬ 
sulted in wrong orbital period determinations (e.g. V404 Cygni; 
Casares, Charles & Naylor 1992). However, in our case, the orbital 
period measured from quiescent spectra is fully consistent with the 
outburst results, supporting CS13 conclusions. Casares (2015) has 
recently presented a correlation between the FWHM of Ha and 
the orbital velocity of the companion star K 2 ~ 0.233(13) FWHM. 
For the case of J1357 this results in a donor star velocity of K 2 = 
967 ±49kms^*. Note that this value is consistent with the inde¬ 
pendent empirical relation derived for quiescent BHT between the 
outer disc velocity (v^) and K 2 (Orosz et al. 1994; Orosz & Bailyn 
1995), which results in K 2 > 806kms^* (see also Torres et al. 
2015). This, combined with CS13 orbital period (fully consistent 
with our results) yields/(Mi) = ll.OiT.lM© {i.e.MBH > 8 . 3 M 0 
at 90% confidence). If we include the BH orbital velocity es¬ 
timated in CS13 {Ki = 43±2kms^*), the obtained mass ratio 
(q = M 2 /M 1 ~ 0.04) provides more restrictive limits to both Mgn 
andM 2 . We findMg// > 9 . 3 M 0 andM 2 > O. 4 M 0 . Here, we use the 
limit value for the orbital inclination i = 90°. This conservative re¬ 
sult places the system as the most massive LMXB BH together with 
GRS 1915-1-105 ( 1 O. 1 ±O. 6 M 0 , Steeghs etal. 2013). Only CygX- 
1, a high-mass X-ray binary, exceeds the BH mass presented here 
(~ 15M0, Orosz et al. 2011). The donor star mass is consistent 
with the constraint on the spectral type (later than M2V; see Sec. 
3.2), which requires M 2 < 0.4M 0 . 
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4.1 Distance and height over the Galactic plane 

The constraint on the veiling factor, combined with the spectral 
type (mass) of the donor star, imposes a more restrictive lower 
limit on the distance to the system. Photometric values were ob¬ 
tained from acquisition images (one per night) in the /-band. We 
use the relation between /, B and V magnitudes described in 
Fukugita et al. (1996), combined with My and B — V tabulated val¬ 
ues of main-sequence stars (Cox 2000) to obtain absolute magni¬ 
tudes, M/. If we consider the latest spectral type for the donor star 
proposed by CS13 (M5V, My = 11.61), we obtain a conservative 
lower limit to the distance of the system. This assumes a some¬ 
what evolved donor from the Smith & Dhillon (1998) empirical re¬ 
lation, instead of a main sequence star. We measure / = 21.6 ±0.2 
from the two nights when the system is faintest, which combined 
with X > 0.81 results in a conservative constraint to the distance 
of / > 2.29 kpc. If we consider the companion star spectral type 
obtained in the previous section (M2V, My = 9.28), the system is 
required to be further than d > 6.7 kpc, a value which is close to the 
upper limit reported in Shahbaz et al. (2013). 

The above constraint to the distance sets new lower lim¬ 
its to the X-ray luminosity during both outburst and quies¬ 
cence. Armas Padilla et al. (2014b, 2013) obtained Tquies = 8 x 
lO^^ergs^* and Tpeak = lO^^ergs^^ considering a distance of 
d = l.Skpc. Our constraint implies: Lquies > 1.9 x 10^*^ ergs^* and 
Lpeak > 2.33 X lO^^ergs^/ These values are still consistent with 
11357 being a very faint X-ray transient and perhaps the faintest 
stellar mass BH known in quiescence (see Armas Padilla et al. 
2014b for a discussion on the topic). 

The lower limit to the distance, combined with the 
high Galactic latitude (b = 50.004°), places the system at 
z > 1.75 kpc above the Galactic plane. Other BHT that 
might be members of the Galactic thick disc population 
(z > Ikpc, Gilmore* Reid 1983) are: SWIFT J1357.2-0933 
(this work and CS13), MAXI J1659-152 (BH candidate, 
Kuulkers et al. 2013), XTE 11118-1-480 (Uemuraetal. 2000; 
Gelino et al. 2006), XTE J1859-I-226 (Corral-Santana et al. 2011), 
H1705-250 (Remillard et al. 1996, Jonker & Nelemans 2004), 
GS 1354-64 (Casares et al. 2004, 2009) and SWIET 11753.5-0127 
(BH candidate, Zurita et al. 2008); see Corral-Santana et al. (2015). 

4.2 On the orbital inclination 

The He I 5876 A emission line exhibits a particular profile with 
a sharp, deep absorption core, even reaching 0.95 the continuum 
value in some spectra. This has to be treated with caution since 
the Nat 5890 A sky line is placed close to the core of the double- 
peaked He I profile, and therefore a deficient sky subtraction could 
affect the profile. In order to minimize this effect we combined the 
data of the two nights were the object was brighter, because this 
allows for a better correction. The final extracted spectrum (see 
Eig. 4) still exhibits a sharp core profile, reaching ~ 1.05 times the 
continuum level. The residual sky substraction of the nearest and 
equally intense sky line [Ol] 5577 A is reduced to almost noise 
level. This result suggests that the deep absorption core is a real 
feature. 

Deep He I line absorption cores have only been pre¬ 
viously observed in eclipsing, high inclination (i ^ 75°, 
Schoembs & Hartmann 1983) cataclysmic variables in quiescence 
e.g. Z Cha, (Marsh, Home & Shipman 1987), which has i = 81°. 
Obscuring material above the plane of the accretion disc has been 
proposed as the origin of the deep, central absorption observed 
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Figure 4. Spectra centred on the He I (5876 A) emission line. Upper spec¬ 
trum: normalised, averaged 11357 spectmm of the two brightest nights after 
sky conection. Lower spectrum: Averaged spectrum of the sky multiplied 
by a constant factor in order to compare with the upper spectmm. 

in He I and Balmer hydrogen lines (see Rayne & Whelan 1981). 
Therefore, it should not be surprising to find similar behaviour in 
high inclination BHTs. However, only systems with inclinations 
up to i ~ 70° have been discovered and detected in quiescence 
so far (e.g. XTE 11118-1-480, Gelino et al. 2006). None of them 
shows such deep absorption cores. Torres et al. (2015) already no¬ 
ticed depth variations in the core of of Ha profile. We also observe 
variations in the normalized flux of the Ha core in the range 1.3 to 
2.1 times the continuum value. 

It has been suggested by Armas Padilla et al. (2014a) and 
Torres et al. (2015) that the X-ray properties of the source (e.g. ab¬ 
sence of emission/absorption lines) and the low extinction, com¬ 
parable to that expected from interstellar material in the line-of- 
sight, argue against the very high inclination (i > 80°) scenario 
proposed by CS13. Given the lack of eclipsing BHs, the closest 
comparison can be made with high inclination neutron star sys¬ 
tems. These only show absorption lines during thermal, soft states 
(Ponti, Munoz-Darias & Fender 2014; Pontietal., 2015). How¬ 
ever, it should be noted that J1357 never abandoned the hard state 
(see Munoz-Darias et al. 2014 for a direct comparison between 
neutron star and BH states). On the other hand, the two eclipsing 
accretion disc corona sources X1822-371 and 2S 0921-63, show 
extinction values consistent with an interstellar origin (i.e no intrin¬ 
sic extinction) when fitted with standard hard state spectral models 
(lariaetal. 2001; Kallman et al. 2003). This suggests that X-ray 
photons emitted in the direction of the outer disc rim do not reach 
the observer, and only those radiated above/below the rim do it, 
solely interacting with the interstellar material in the line-of-sight. 
Therefore, we still consider i > 80° as the most plausible value. 
The detection of He I line cores and the already large BH mass (for 
i = 90°) presented here further support this conclusion. 


5 CONCLUSIONS 

We have used 10.4-m GTC optical spectra of Swift J1357.2-0933 
to constrain the accretion related contribution to the optical emis¬ 
sion in quiescence. This yields more restrictive lower limits to 
the distance (d > 2.29kpc) and height over the Galactic plane 
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(z > 1.75kpc), placing the system in the Galactic thick disc. We 
detect variability in the Ha profile modulated with a period P = 
0.11 ±0.04d, confirming the period detected in outburst. Using 
the recently presented FWHM — K 2 correlation, new constraints 
to the fundamental parameters are derived. In particular, we obtain 
Mbh > 9.3M0. We also favour q ~ 0.04 and M 2 O.4M0. This 
indicates that Swift J1357.2-0933 harbours one of the most massive 
BHs known in our Galaxy. 
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